THE PROFILE OF BUBBLING SOLUTIONS 
OF A CLASS OF FOURTH ORDER 
GEOMETRIC EQUATIONS ON 4-MANIFOLDS 

GILBERT WEINSTEIN AND LEI ZHANG 

Abstract. We study a class of fourth order geometric equations defined 
on a 4-dimensionaI compact Riemannian manifold which includes the Q- 
curvature equation. We obtain sharp estimates on the difference near 
the blow-up points between a bubbling sequence of solutions and the 
standard bubble. 



1. Introduction 

Let (Af , g) be a compact Riemannian manifold. The conformal class of g 
consists of all metrics g = e'^^g for any smooth function u. A central theme 
in conformal geometry is the study of properties that are common to all 
metrics in the same conformal class, and the understanding and classification 
of all the conformal classes. For this purpose it is often useful to be able to 
single out a unique representative in each conformal class by imposing some 
geometric condition. This usually leads to a conformally covariant geometric 
equation for the conformal factor e^". Such equations have attracted much 
interest in the literature in the past half-century. 

In dimension 2, the natural condition to impose is constant Gauss cur- 
vature. The Poincare Uniformization Theorem states that this is always 
possible: every compact Riemannian surface is conformal to one with con- 
stant Gauss curvature. The conformally covariant operator in this case is 
the Laplacian-Beltrami operator, given in local coordinates by: 

and the equation for constant curvature is: 

-AgU + Kg = HgC^"^ 

where Kg is the Gauss curvature of g, and Kg is constant. Using the Gauss- 
Bonnet theorem, we see that the sign of Kg is determined by xm the Euler 
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characteristic of M: 

27rxM = / KgdAg = Kg / e^'^dAg, 
Jai Jai 

where dAg is the area element of g. Although this result was originally 
proved by Poincare using non-PDE methods, there is now a PDE proof, 
see [TT| and |37j . Furthermore, the operator Ag is conformally convariant 
= e-^-Ag. 

For compact Riemannian manifolds of dimension n > 3, a natural gener- 
alization is to impose constant scalar curvature. This leads to the Yamabe 
problem: given a compact Riemannian manifold (M, g) of dimension at 
least 3, find a metric conformal to g with constant scalar curvature. This 
was also eventually answered in the affirmative, see |39l [34l [21 [32]. The 
corresponding operator is now the conformal Laplacian Lg = Ag — CnRg, 
Cn = {n — 2)/4(n — 1), and the equation for constant scalar curvature is the 
Yamabe equation: 

LgCj) = e<^(-+2)/(n-2) 

where (/)^/("~2) _ g2n^ ^^^^ ^ jg or —1. The operator Lg also has a 

conformal covariant property: 

where g = g. 

In 4-d, another problem analogous to the 2-d case arises from imposing 
the condition of constant Q-curvature: 

Qg = -l^{AgRg-Rl+?,\mCg\''). 

The natural question is the same: given a A-d compact Riemannian manifold 
(M,g), is there a metric g = e^'^g in the conformal class of g with constant 
Q-curvature? The Q-curvature of the metric g is given by: 

PgU + 2Qg = 2Qge^^, 

where Pg is the Paneitz operator: 

Pgu = Alu + divg (^(^^Rg9 - 2Ric3^ Vu 

Integrating with respect to the volume element dVg, it is easy to see that 
the quantity: 

kp = Qg dVg 

Jai 

is a conformal invariant, i.e., it is constant in the conformal class of g. 
Furthermore, we also have a Gauss-Bonnet formula: 



M 



Q9 + I \Wg\^ ) = 4vr'xM, 



ASYMPTOTIC EXPANSION FOR BUBBLING SOLUTIONS 



3 



where Wg is the Weyl tensor of g given in local coordinates as: 

2 2 

Wijki = Rijki - :^^^i9i[kRi]j - 9j[kRi]i) + (^_;^)(^_2)-^^«['=^'li- 

We note that Wg is pointwise conformally invariant Wg = Wg, and the 
operator Pg is conformally covariant: 

(1-1) Pgf = e-^^Pgf. 

We use ([P]) to denote the assumption: 
{P) Ker(Pg) = {constants}. 

We remark that is often satisfied. For example, Gursky in [T7] proved 
that if (M, g) has non- negative Yamabe invariant Yg>Q, and satisfies kp > 
0, then ([P]) holds and Pg > 0; see also [TH] where the assumptions are 
weakened to Yg > and kp + /6 > 0. 

Chang and Yang proved in [9j that \ikp < Stt^, Pg > and ([P|) holds, then 
there is a conformal metric g whose Q-curvature is constant. In [13j, Djadli 
and Malchiodi extended this existence result assuming only that ([P|) holds 
and A;p / Svr^A^ for any positive integer N . An essential ingredient in this 
existence result is an a priori bound: if kp ^ Stt'^N for any positive integer 
A^, then any sequence of solutions of the prescribed Q-curvature equation is 
uniformly bounded. In fact, this a priori estimate can be extended to the 
following more general equation in the same class: 

(1.2) Pgu + 2b = 2/ie^", 

where 6 is a smooth function. Note that if 6 = Qg, then h is the Q-curvature 
of the conformal metric e'^'^g. Assuming ^ Hq, > cq > 0, and ^ bo, 
Druet and Robert in [151 showed that any sequence of solutions {uk} of (|1.2p 
with h = hj^ and b = b^ is uniformly bounded, provided J^j 6o 7^ Stt'^N, see 
also Malchiodi [25], . 

However, bubbling can occur when J^^ bo dVg = Sir'^N for some positive 
integer A^. A precise understanding of this bubbling phenomenon is required 
if progress is to be made on the existence problem. The study of the blow- 
up profile and other blow-up phenomena for the Paneitz operator and other 
4-th order elliptic equations has attracted much interest recently; see for 
example Pia[8l[l2l[l6l[ia[22l|26l[27l[28l[30l[Ml[M^ 

Let {uk} be a sequence of solutions of (|1.2p with h = h^, and b = bk- 
We say that this is a bubbling sequence if suplu^l oo. In [15], Druet 
and Robert studied bubbling sequences of solutions of (|1.2p and obtained 
some asymptotic estimates on the behavior near the blow-up points. We 
will throughout make the following assumptions on the coefficients bk and 
hk- 

{b, h) \\bk - boWc^iM) 0, \\hk - /io||c2(Af) ^0, hk> cq. 

It follows immediately that ||6a:||ci(m) ^ Co, and ||/ifc||c2(M) ^ C'o for some 
constant Co independent of k. We let G denote the Green's function for the 
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Paneitz operator: 

(1-3) f{0-f9=[ G{^,r,)Pgf{v)dVg{v), [ Gi^,ri)dVg{v)=0, 

JM JM 

where fg = \o\g{M)-^ f^j f dVy is the mean value of /. The asymptotics of 
this Green's function are studied in the Appendix. Now, for fc = 0, . . . , let 

(1-4) 0fe(O = 2 / G{^,ri)hk{'n)dVg{ri). 

JM 

Since {life} is a bubbling sequence, it follows immediately that Jjy^bodVg = 
SNtt'^ for some positive integer N. Druet and Robert proved that passing 
to a subsequence, there is a finite set S = {pi, ..,pn} such that: 

TV 

Uk-Uk^ 167r2 ^ G{pi , •) - 00 in C^M \ S) , 

i 

Let P be the regular part of the Green's function: 

(1.5) G(e,7?) = --^Xir) logdg{C,ri)+/3{C,r)). 

Here x is a cut off function supported in a small neighborhood of ^, and 
r = dg{^, T]). They also proved that for z = 1, . . . , A?^: 

647r2V2/3(K,K) + 647r2 ^ ViG{pi,pj) - 4V0o(k) = '^JT^^ 

where h is the limit of hj^ as k ^ oo, and Vi, V2 denote the derivatives 
with respect to the first and second variables respectively. In this article we 
will continue this line of investigation and derive more precise asymptotic 
estimates for the behavior of such solutions. We define the standard bubble 
at p: 



HdgiMY 

4\/3e 



We will also adopt the following notation. For k large enough, there are 
N points {qik} such that qik Pi and Uk{qik) 00. Let Hik = hkiqik), 
£ik = e-"'=fe'=), and Uik = Uq.^,eik,Hik- 

Theorem. Let {u^} be a bubbling sequence of solutions on M. Then passing 
to a subsequence, there is a constant S > such that for any fixed r G (0, 1), 
there exists a constant Ci = Ci{N,g,co,Co,T) such that: 

(1.6) \ukiO - UikiOl < Cidg{qik,0^, 
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in B{qik, S), and such that for i = 1, . . . , N we have: 



(1.7) 64TT^V2/3{qik, Qik) + 647r2 ^ ViG{qik, q,k) 



4V(/.fe(Qifc) + 



hk{qik) 



< Cie 



Our approach is motivated by Lin and Wei's work [25 from which one 
can easily derive an 0(1) bound in (II. 6p (i.e. \uk — Uik\ < C near pi) pro- 
vided (M, g) is locally conformally flat; see also [38j which uses a completely 
different approach. Our result removes the hypothesis of local conformal 
flatness and also improves the estimate near the blow-up points. We hope 
our approach can be fine-tuned to yield better yet estimates as required to 
handle the existence question posed above. 

A major difficulty when trying to prove a priori estimates for solutions of 
fourth order elliptic equations is the lack of a maximum principle. In order to 
remedy this, Lin and Wei devised a strategy based on the Pohozaev identity. 
We adapt this approach to the case in which the manifold is not necessarily 
locally conformally flat, making use of conformal normal coordinates. These 
are normal coordinates for a metric g in the conformal class of g for which 
det{g) = 1. The existence of such a metric is proved in [6]. Although, we 
used this result for the sake of simplicity, our proof only relies on the weaker 
concept already introduced in the solution of the Yamabe problem where 
one only requires det(^) = 1 to hold to high enough order in the distance 
from the center of the ball under consideration, see |21j . 

We now briefly sketch the outline of the paper and the proof of our The- 
orem. 

In Section [21 we prove the 0(1) estimate. We use the Green's represen- 
tation formula, together with rough estimates from [15j, to write long range 
asymptotic formulas for the rescaled solution and its derivatives in terms 
of the concentration of energy near the singular point, i.e. within a care- 
fully chosen radius Ik = —Eklogek, where Sk is related to the maximum of 
Uk, see (12. 9p . These are then substituted into an asymptotic Pohozaev iden- 
tity, and after estimating the higher order terms, we obtain an asymptotic 
formula for the energy Uk ~ IGvr^. When substituted back into the asymp- 
totic formula for v^, this yields a long range 0(1) bound. Finally, we use 
standard estimates in the interior, and then these long range and interior 
estimates in conjunction with the maximum principle in the mid-range. 

In Section [31 we prove (jl.Gh by contradiction. We divide the argument into 
two cases, depending upon whether an appropriately weighted supremum 
runs off to infinity or remains in a bounded region along a subsequence. In 
the first case, we use the Green's representation formula and a comparison 
between the geometric and Euclidean distances to reach a contradiction. 
In the second case, we show that the difference between the appropriately 
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rescaled solutions and the standard bubble converges to a solution of the 
linearized equation which we can then show vanishes thanks to a lemma of 
Lin and Wei from [24J, again leading to a contradiction. 

Finally, in Section lU we use our estimate ()1.6p in the Pohozaev Identity 

— 1/2 

over a ball of radius to obtain a Euclidean version of the vanishing 

rate. We then translate this result into the original metric g and prove (jl.7p . 

The Appendix deals with delicate estimates for the Green's function, an 
asymptotic comparison between the geodesic distance and the Euclidean 
distance in conformal normal coordinates, some well known curvature and 
metric derivatives computations in conformal normal coordinates, and a 
proof of the asymptotic Pohozaev identity. 



2. The 0(1) estimate 
In this section we derive the 0(1) estimate, i.e., we show that 

WkiO - UikiOl < C, for ^eB{q,k,S). 

This estimate has been established by Lin- Wei [23] for locally conformally 
flat manifolds; see also |38j for a completely different proof. Furthermore, 
we remove the assumption that {M,g) is locally conformally flat. 

Our first step is to rescale the solutions, and use the Green's represen- 
tation formula (jl.3p to derive the long range asymptotic formulas (|2.1ip - 
(I2T5D . 

In [15], Druet and Robert prove that the singular set S consists of only 
finitely many points {pi, . . . ,pjv} and these are separated uniformly in k by 
a positive distance. Without loss of generality we will focus in this section 
on pi, and to simplify the notation, we will omit the subscript 1, so that we 
now consider a sequence of points qk & M where Uk has a local maximum 
UkiQk) — > oo and — > p as ^ oo. 

According to [6], we can find function Wk defined on M, such that in a 
neighborhood B{qk,5i) of qk, Si > 0, we have det(^fc) = 1 in the normal 
coordinates of the conformal metric gk = e^'^'^g. We refer to these coordi- 
nates as conformal normal coordinates. We point out that det(^fc) ~ 1 to 
high enough order would be sufficient for our purpose, but we use Cao's 
result since it simplifies the proof. We also choose 5i small enough so that 
5i < inj(M)/10 and 6i < d/10 where d is the minimum distance between any 
two points in the singular set S. Using the conformal covariance property 
of Pg (jl.ip . we obtain that the function Uk = — Wk satisfies 

where 2hk = Pg^Wk + 26^6^^"''=. We remark that if hk = Qg then bk = Qg^- 
We also note that WkiC) = 0{dg{^,qk)'^) in a neighborhood of qk, hence all 
the terms coming from Wk can be absorbed on the right-hand side of (|1.6p . 
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We have the fohowing estimates, also proved in [15j : 
(2.1) 



\D'uk{i)\<C{dg,{i,qk))-', j = 1,2,3, 

where \D^Uk{£^)\ = |Z)'^{tfc(^)| and the sum is over all multi-indices J of 
order j, and Ek = e^"'=*^'^'=\ We now rescale the solutions Uk, using a blow- 
up of the neighborhood of the point qk- Define the map (pk'- -6(0, (^le^^) 
B{qk,5i) hy tpk'- y ^ SkV, where on the right-hand side we are using confor- 
mal normal coordinates on B{qk,6i). We use the notation f = = f o ip 
to denote the pull-back of a function / defined on B{qk,6i), and we let 
gk = £k'^'P*gk be the blow-up metric, i.e., a rescaling of the pull-back met- 
ric. We define: 

Vk = Uk + logefe, 
and note that Vk{0) = 0. It follows from (jl.ip that Vk satisfies: 

(2.2) P^^Vk + 26% = 2^6^'" , in B{0, 6ie^^). 
The estimates (12. ip now read: 

(2.3) |t;fc(y)| < (-2 + /X) log(l + |y|) + C7(/i), |y| < ^i^^^ 

(2.4) \Divk{y)\<C{l + \y\)-^, j = 1,2,3. 

where /i G (0, 1). 

Let Gk be the Green's function for P^^ . Then, we have: 

Uk{0=^k + 2[ Gk{i,'n)hk{ri)e^''^^^Uv^,^{r,)-2 I Gk{tv)bkin)dVg,{v)- 



M JM 

where Uk is the mean value of Uk- Decompose Gk into a principal part and 
a regular part as follows: 

Gk{i,r]) = -^x{r) logd9,(e,^)+/3(e,r/) =F(e,r/)+/3(C,r?) 

where x = 1 on B{qk,6i), x = on M \ B{qk,25i). We have 

(2.5) Uk{0=^k + 2f H{^,r,)hk{v)e^'''^^^UVg^{rj) + 4>k{0 

JM 

where 
(2.6) 

4(0 = 2/ P{tv)hk{v)e^'"=^'^UVg^{v)-2 f Gk{^,v)h{v) dVg^{v). 

JM JM 

Note that since det{gk) = 1 in B{qk,Si), we have dVg^{7]) = drj in B{qk,6i). 
Taking the difference of ()2.5p evaluated at ^ and pk, we get: 

(2.7) n,(^)-^,(p,) = ^^^log(i^^) Xir) hkiv) e'^'^^'^UVgM 

+ 4(0 - 4(9fc)- 
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Here we have used the fact that smce the coordinates are normal dg^ (t], qk) = 
\r] — qk\. Thanks to the cut-off function Xi '^^ can now replace the integral 
over M by an integral over B(qk,26i), and after rescaling, we now obtain: 



(2.8) Vkiy) = / log ( , xiekr) hk{ekz) e^^'^^^) dz 

+ ^k{£ky) - (i>k{o)- 



Let 

and define 
(2.9) 



Zfc = -efelogefc, Lfc = -logefc, 



By (j2.3p and the fact that det{gk) = 1, one sees easily that: 
(2.10) ak = 2 [ hkiv) e^"^^^^^ di^ + 0{Lf). 

As in |24j, the representation formula (j2.8p implies the following long 
range asymptotic formulas for Vk and its derivatives: 



(2.11) 


Vk{y) 


= -01og|y|+O(l), 


Lk < \y\ < he^^ 


(2.12) 


drVk{y) 


= -^Lk' + OiL-% 


\y\ = Lk 


(2.13) 


dr{rdrVk{y)) 


= 0{Lf), \y\=Lk 




(2.14) 


^Vkiy) 


= -^Lf + OiLf), 


\y\ = Lk 


(2.15) 


drAvk{y) 




\y\ = Lk- 



Note that while the asymptotic formula is required in the whole range Lk < 
\y\ < ^i^k"^ foi' ^fc) it is only required on \y\ = Lk for the derivatives. 

Since the proof of these estimates is similar to the one in [2l|, we will 
only briefly sketch the argument pointing out the main differences, a major 
one being the difference between the Euclidean and Riemannian distance. 
It follows from (j2.8p that for \y\ > Lk'- 

(2.16) vk{y) = --^ [ logd^g^{y,z)hk{ekz)e^^^^'Uz + 0{l). 

Moreover for any multi-index J of order j = 1,2,3, and \y\ = — loge^, we 

have: 

(2.17) 

D'vkiy) = - A / D^{logd^g^{y,z))hk{ekz)e^'"'^'Uz + 0{ei). 
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We now divide the domain of integration in ()2.16p and ()2.17l) into three 
subsets B{0, Sis'i^^) = Jli U U ils, where: 

f]i = {|z| < |y|/2}, n2 = {\z-y\<\y\/2}, = B{0,6ie^')\{QiUn2). 

We also use the fohowing approximations of the distance and its deriva- 
tives by their Euchdean counterpart^: 

(2.18) logd§,(y,z)-log|y-z| = 0(1), z € B{0, 8^^5i/2) 

(2.19) \D^{\ogdg,{y,z)-log\y-z\)\<Cel\yf~^, z e fti. 

Over r22Ur23, the integral (I2.16P can be estimated simply by using ()2.3p and 
the approximation ()2.18p leading to: 

/ logd^^{y,z)hk{z)e^^>^^''Uz = 0{\y\-^+f'^), /xi > small. 

In order to capture the asymptotics of the integral ()2.16p over f2i, we again 
use the approximation ()2.18p to reduce the calculation to the Euclidean case, 
so that (pflT]) then follows with the help of ([2^0]) . 

Similarly the estimate of (I2.17P over U ils can be obtained from the 
bounds: 

\Dnogd^g^{y,z)\<C\y-z\-^, j = 1,2,3, zGl^aU^a, 
leading to: 

/ \D^ {log dg,{y,z))\ h{z)e''^'^^^Uz = 0{\yr'-^+n, 

while the estimation of these integrals over Qi requires the more precise 
approximation ()2.19p . 

We now use the long range estimates (j2.1ip - (|2.15p in the following Po- 
hozaev identitjfl over the ball 17 = B{qk, 1^): 

Jn ^ JdVL \^ 

+ ^ (^Agiidig'^dju + CAgudimg^'djU + C^gud^m9''di,u - 2bediU 

+ 2 / (R,j,i{0)C^C^idjud^u + O{r^)\Du 
Jdn ^ 

2Rij^i{0){^^djudiU + CC^djudi^u) + 0{r^)\Du\^ + 0{r^)\D^u\ 



n 



Here, we used the conformal normal coordinates on this ball, we denoted 
r = 1^1, and denoted the unit normal to the boundary by Ui. Furthermore, 



^ These approximations are proved in the Appendix 
^The proof of this identity can be found in the Appendix. 
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to simplify the notation, we suppressed the sequence index k, and since 
det(g) = 1, we omitted dVg = d^. Finally, we remark that we chose to write 
this identity in terms of Uk rather than Vk to avoid an even longer formula. 
It is easy to translate the long range estimates (]2.1ip - (l2.15p to iik from the 
fact that Vkiv) = Uk{eky) + logefc. 

We denote the integral on the left hand side of this identity by Iq) and 
the four integrals on the right-hand side by Ii, I2, Is and I4 respectively. 
By (j2.3p . we obtain: 



2 

hence it follows from (j2.10p that: 

(2.20) Io = ak + 0{Lf) 

By the expansions (jC.lip and ()C.12p of the derivatives of the metric 

and (j2.4p we get: 

(2.21) 

\I2\<C [ el\D\k\ \Dvk\ |yp + el\y\^\D\k\^ + et\Dvk\ \y\ = O(efc) 

and similarly, using ()C.14p . we see that: 

(2.22) II3I + I/4I = O(efc). 

It remains to compute Ii. First, using (|2.3p . we can estimate the first 
term in Ii: 

i / e^^hke^'''= = 0{Lf). 
^ Jan 

Using this bound, and using the expansions (jC.OP and (jC.lOP in the remain- 
ing terms, we can now reduce Ii to: 

Ii = y (^lkdu{l^Uk)duUk + AukduUk + ViCl^UkdimUk - ]^lk{Aukf^ 

+ 0{ek) 

= [ ( -Lkdu{Avk)duVk + duiy ■ Vvk)Avk - ^Lk{Avkf] + 0{ek). 
Using (I2.1ip - (12.15p in the above, we get: 

(2.23) ^i = I^ + ^(^^"')- 
Combining (I2:2TD . 0:22]) and (l2:23D . we get: 

»k + 0{Lk') = -^,+0{L,^). 

which implies 

(2.24) afc = 16^2^ 0(L^^). 
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When substituting this into ()2.1ip . we obtain: 

Vkiy) + 2log\y\=0{l), \y\>Lk. 

The argument in the region \y\ < Lk follows the one in [_2.4J closely, hence 
we again only sketch the proof. Without loss of generality, we assume that 
hkilk) fl Let U = f/o,i,i be the standard bubble in M^: 

It is easy to check that U satisfies A'^U = 2e^^ and it is well known that 
Vk ^ U in Cj'^^(M^), see for example [TB]. Thus, for any fixed A and all k 
sufficiently large, we have: 

\vkiy)-U{y)\<l, for|y|<A 

Subtracting (j2.17p from its Euclidean counterpart, using (j2.19p to compare 
d^{y, z) and \y—z\ as well as their respective derivatives, and also using (j2.24p 
to compare the leading terms, we obtain: 

\Avk{y) - AU{y)\ < C\y\'\ for A<\y\< Lk- 

Now, letting T{y) = C(l + and choosing C large enough, we can 

guarantee that AT < —\/\vk — Af7| whence from the maximum principle 
\vk{y) - U{y)\ < C(l + \y\~^), on ^ < \y\ < Lk- Substituting ^ = e^y, and 
using the definition of Vk, we obtain the version of (jl.6p with 0(1) on the 
right-hand side, i.e., with r = 0. 

3. A Sharper Estimate 

The main purpose of this section is to establish (jl.6p . An important tool 
we use is the following lemma, due to Lin- Wei [24J: 

Lemma 3.1. Let U{y) = -log(l + |yp/4\/3) be defined on M^. Then U 
satisfies 

A^U = 2e^^, C/(0) = maxC/ = 0. 
Furthermore, any solution of the linearized problem: 

A^cP = 8e^^0, |0(y)| < C(l + r G (0, 1), 

is given by (p = Yl^=o ^i"^] where 

^ 1 - |y|V4V3 
l + |y|2/4^/3 

l + |y|2/4\/3' ' ' 

Remark 3.1. One immediate consequence of Lemma \3.1\ is that if in addi- 
tion, (j) satisfies (piO) = 0, and V0(O) = 0, then (p = 0. 

■^Otherwise, we can add a constant to fife. 
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Let pk = /iA:(0)^/^/4\/3 and consider the solution f/^ (y) = —\og{l+ pk\y\^) 
of the equation 

(3.1) A^C/fc = 2/ifc(0)e^^^ t/fc(0)=0, |VC/fc(0)| = 0, 

on M^. Letting Wk = Vk — Uk, then by the result of Section O we already 
know that \wk\ < C in B{{),5ie~^^). Our goal in this section is to prove: 

(3.2) \wk{y)\<Cel\y\\ \y\<hek^. 

for any < r < 1, which implies (|1.6I) . Equivalently, if we let 

\wk{y)\ 



max 



el{l + \y\ 



then it suffices to show that is bounded on Q.^ = B{{),5ie^^). Suppose 
that Afc — > oo, and let yt € ^k be the point where A^ attains its maximum. 
Now, either: (i) yk — > oo; or (ii) \yk\ remains bounded at least along a sub- 
sequence, and hence a further subsequence, which without loss of generality 
we will assume is yk itself, converges to y* . We will show that in both cases 
a contradiction follows. 
Define: 

- / X wk{y) 
^^^y> = A ^rn _L I — ivT- 

Afc4(i + \yk\) 

By the definition of A^, we have 



(3.3) \wk{y)\ < 



1 + \yk\ 



and Wk{yk) = ±1- 

Assume first that yk ^ oo. Since \wk{y)\ < C and A^ — > c«, we clearly 
have yk = o(l)e^^. /,From the fundamental solution for A^, it is straight- 
forward to get: 

(3.4) ' I ' 

= / log hk{0) e4^'=(^) dz + 0(4). 

Similarly, using the fundamental solution for P^^ , we find: 

My) = A / log T^T^ ^kiskz) e^^''^'^ dz + 0{ek\y\) 

^^■^^ 1 r 

= T^ / hkiekz)e^'^>'^'Uz + 0{ek\y\), 
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see (I2.8p . Note that for the second equahty, we used (jB.Qh as well as the 
decay rate of Vk- Finally, we estimate the source term: 



= \hu{euz) {e^^^ - e^^'=) + - /i,(0)) e^^^' | 

<C{l + \z\r''\wk{z)\+0{ek){l + \z\)-' 

< CelAkil + + Oiek){l + \z\r' . 

Substituting this in (|3.5p and combining with p.4p in the definition of Wk, 
we obtain: 

WkiVk) = / log ■ 



fife V + 

At(l + I»l)' y * 

Since oo, it is now easy to see that the right hand side is o(l), which 

contradicts Wk = ±1. 

We now turn to the second case and assume without loss of generality that 
Uk converges to y*. We will show that along a subsequence Wk converges. 
This will be accomplished by estimating Pjj^iUk — Vk)- We start with: 

PgUk = Apk + 4 divg (^(^^R{eky)gij{y) - 2Rij{eky)^ dU^ 

where we have suppressed the subscript k on the metric and curvature com- 
ponents, and where R, Rij are the scalar and the Ricci curvatures of g. In 
conformal normal coordinates, we have: 

(3.6) Apk = A'^Uk in B{0, die^^^). 

Furthermore: 

= dn.r' (^lRi^ky)h - 2R^jieky)^ g'^diUk 
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Since = = 0(4lyl')> ^i^kV) = 0{el\y\''), and 

Riji^kV) = 0{sk\y\), we easily get Ai = O(e^), and A2 = 0{sl). Further- 
more, since in addition Rij^i{0) = — ^-Rj(O) = 0, we also have = 0(e|), 
and A4 = 0{e1). Substituting this into the above equation and multiplying 
by el, we get: 

(3.7) eidiv^ (^(^lR^eky)gi,iy) - 2Rij{eky)^ dU^ = 0{el). 



Combining (j3.6p and (j3.7p . we find: 

(3.8) PgUk = 2/ifc(0) e^^*^ + 0(4), \y\ < diS^K 
Combining (13. 8p with (12. 2p . we obtain: 

(3.9) P^«;fc = 8/ifc(efey) e^^" wu + 0(efc)(l + |y|)~^ + 0(4), |y| < (^^le^ \ 

where S,^ is given by: e^^*" = Jq*^ g4tofc+4(i-t)c/fc Finally, this leads to the 
following equation for w^: 



(3.10) Pgta, = 8/^fe(6,;,)e^g''^fc + ^^'^ 1^ 'it 'Ah^l IV 

Afc(l + |yfc|)^ Afc(l + |yfc|)^ 

Since — > y*, a subsequence of converges to w in C^(IR^). We will 
assume without loss of generality, as in Section [2l that /^^(O) ^ 1. It follows 
that the limit w satisfies: 

\w{y)\<C{l + \y\r, 
w{0) = \Vw{0)\ = 0. 

By the remark following Lemma [3.H we conclude that w = 0, which contra- 
dicts w{y*) = ±1. This concludes the proof of (j3.2p . 

In Section m we will also need estimates on the derivatives D^Wk{y) for 
\y\ < ^fc ^1^1/2, j = 1,2,3. By combining (|3.9p and (j3.2p . we have: 

Pgw,{y) = 0(4)(1 + \y\)-'^' + 0{et), \y\ < ^-^e^\ 

Fix y, and let r = \y\ and fkiz) = Wk{rz) for 1/2 < \z\ < 2. Then fk{z) 
satisfies: 

P^fkiz) = 0(eD(l + r)-^^^ + Oiey), B2 \ B,/^, 
fk{z)=0{elr-), B2\By2, 

where g is the rescaled metric r~'^ip^g and ip: z ^ rz. Standard elliptic 
theory for fourth order equations [5j yields: 

D^fk{z) = 0{elr-), \z\ = 1, j = 1,2,3. 

Hence, we conclude: 

(3.11) D^wkiy) = 0{el\yrn, 3 = 1,2,3, \y\ < e7^6i/2. 
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4. The Vanishing Rate 

The purpose of this section is to complete the proof of our Theorem by 
proving (II. Th . We first prove a Euchdean version in Subsection l4.1|, and then 
translate this result to the original metric in Subsection 14. 2[ 



4.1. A Euclidean Version. The goal of this subsection is to prove (I4.14|) . 
This is accomplished in three steps: 

(i) In the first step, we derive an asymptotic expansion of the Pohozaev 
identity; see (j4.3p . 

(ii) In the second step, we express this identity in terms of the back- 
ground Euclidean metric; see (j4.12p . 

(iii) In the last step, we complete the proof of ()4.14p . 



— 1/2 

4.1.1. Step 1. In this subsection we let Ek = i?(0, e^, ), and we derive 
an asymptotic Pohozaev identity for Vk on E^. We multiply (12. 2p by daVk, 
a = 1, ..,4, integrate with respect to the Euclidean volume element dy, and 



estimate each of the resulting terms. First by the 0(1) estimate and { b,h ): 
(4.1) / 2hkieky)e^^''(y^daVk{y)dy 



^ -'Ek ^ JdE. 



= -'-^dMo) [ e^^'=(^)+o(4: 

J El; 

Next, integrating by parts, we have: 

(4.2) / AlvkdaVk 
Je, 



k 



where we used: 

djig'^diaVk) = da{AgVk) - diag'^djVk - dag'^dijVk 

We now estimate the two integrals over Ek in (j4.2p above. Using diag^-' 
0{el\y\), which is implied by (jClil) . and (I23D, we find: 

A-gVk diag'^ dj vk = 0{el). 



E, 



k 
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Furthermore, for any < o" < 1, the second integral over Ek in ()4.2p can be 
estimated as follows: 



(A^vk) dap dijVk 
= 4j^ (AUk + AijWk) (^-'^Ri(a«^)jiO)y"' + 0{ek\y\'^)J dijVk 

where we used (j3.1ip , and the following estimate implied by (jC.lip : 

as well as the antisymmetry of the curvature tensor, and Rij{0) = 0. Here, 
we use the customary round brackets notation to denote the symmetric part. 
We will choose a E (0, 1) at the end of the argument of Section (j4.ip . Next, 
since R{eky) = 0(e^|yp), and Rij{eky) = 0{ek\y\) and ([23D, we have: 

elj^ dm [r' {^R{eky)9i3 - IRiA^m)^ dlVk (}^^ daVk 

= el I 5™ ( \R{eky) mj - 2Rij{eky)] diVk daVk Vm 
-elj^ h - 2Rij{eky)^ OiVkg'^ OamVk = 0{el). 

Finally, we estimate: 

el [ 2bkdaVk = 0{el^^). 

JEk 

Combining all the terms, we arrive at the following Pohozaev identity: 
(4.3) '^dahkiO) [ e^-^dy + 0(4'+'^)/') 

^ JEk 

= j (^-g'^ dj{AjVk) daVk Vi + A-gVk diaVk - ]^{A"gVkf ■ 

4.1.2. Step 2. In this second step, we rewrite ()4.3p in terms of the Euclidean 
Avk rather than A^Vk- We begin by substituting: 

iPiy) - n = -\elRiimj{0) ymyiY^ + 0{el\y\^) = 0{el\y\^). 

"J \y\ 
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into (Ol) to get: 



(4.4) '-^dahml e^^^dy + 0(4='+'^)/') 



dEk \ ^ 

A straightforward computation leads to: 

(4.5) di{A^gVk) - di{Avk) 

= din.g'^^ divk + diivk + dig"^^ dmiVk + (5""' - ^mi) diM- 

1/2 

In view of ()C.14p and ()C.12p . we have, for \y\ = : 

(4.6) dimr' diVk = Oiel), dmr' duvk = 0{el), 

(4.7) d,r' dmiVk = d,r' {d^iUk + 0{elr'^-^)) = 0(ef +'^)/'), 
where to derive ()4.7p . we also used the following consequence of (jC.lip : 

di^' = -\Rm(ia)melya + 0{el\y\''), 

as well as the anti-symmetry of the curvature tensor and Rij{0) = 0. Next, 
on lyl = Ef, ^^'^ , we have 

(4.8) ir'' -S'^')^^mlVk 

lelRmabliO) Va Vb + 0(4^')) {dimlUk + 0{elr''-'')) 

= - ^4^^aM (0) [/, + O (ef ) = O (5(^+^)/^ ) , 

where we have used the following expansion, valid for any radial function 
f{r): 

5-/(0 = (r (0 - ^ + ^) ^ 

+ [ f"{r) - Vm + yi Sim)r^ - 2yi ym Vi 



+ {f"{r)-f'{r))- 



as well as the anti-symmetry of Rabcd and Rij{0) = 0. It now follows 
from (113]), (1121) and (gSD that the following holds: 

(4.9) - / d^{A^gVk) daVk yi = - I di{Avk) daVk Vi + 0(4'+"^^'). 

JdEi, JdEu 
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Next, on \y\ = e^^^^, we have: 

^gVk = dif^ djVk + g'^ dijVk = 0{el)r + Avk + {g'^ - 5ij) dijVk 
= 0{el) + Avk - (^^RiabjiO) Va Vb + 0(4^^')) [dijUk + 0{elr^-^)) 

from which it follows: 

(4.10) / AgVkdiaVkVi= [ AvkdiaVki^idS + 0{e^,^^''^^^). 

JdEk JdEk 

Similarly: 

(4.11) -U {A^vkfva = -\f {Avkfva + 0{ef^''^"). 

Substituting (fOj) . (fiJOj) and (fill]) into the Pohozaev Identity (fO|) . we 
conclude: 

(4.12) '-^dMO) [ e^^'=dy + 0(ef+'^)/') 

= / ( -di{Avk)daVkiyi + AvkdiaVki^i - l:{Avk)'^ Ua] dS. 

JdEk V ^ J 



4.1.3. Step 3. In this subsection, we first aim to replace Vk by Uk in the Po- 
hozaev identity (j4.12p , after which many of the terms will simplify thanks to 
the radial symmetry of Uk , leading to the Euclidean version of the vanishing 
rate (|4.14|) . Recall the definition of (pk (|2.6|) . from which we have: 

D^4>ki^ky) = 0{ei), j = 1,2,3. 

For \y\ = z^^"^ , we cut -6(0, (5ie^^) into three subdomains -6(0, (5ie^^) = 
r2iUr22UJ73 as in Section [2l page[9l and use the representation ()2.8p . Using 
standard estimates over U and (jB.ip over l^i, we find: 

daVk{v) 

= ~T2 I da[\ogdg{z,y))hk{z)e^''^^'''^+ekdaM^)+0{el\y\) + 0{\y\-'') 
= ^ 1^ hkiskz) e^^'=(^) + Ek daW^) + 0{ef), 
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where we have omitted the standard volume element dz. Similarly, 

2vr^ JvLi \y- 



d.My) = /_ ' 1:^^:14 " ^ h,{e,z) e-^ + 0(4). 



\y-A^ 



7^1 \y - zr 

We also have 

= /ifc(0)e^^'=(^) + O(4)(l + |z|)-8+^ 
Substituting this in the derivatives of we have: 

/^i \y-z 



47r^ /o, y — 2r 



^vk{y) = - A / ^A.(0) e^f^'^^^) + 0(e 



o-+l^ 



a.(At;,(y)) = / hm e'^^ + 0(4+'^^)- 

We perform a similar computation for and take the difference, leading 

— 1/2 

to the following estimates for |y| = : 



ai(At;,(y)) = d,{AUk){y) + 0(e^+=^/'), 
Ar;fc(y) = A[/fc(y) + 0(e^+ 
dijVkiy) = 9^jUk{y) + 0{ep 



k 

o-+l\ 



Substituting these estimates into the Pohozaev Identity (j4.12p . we obtain: 
(4.13) '-^dahUO) [ e^^^+0(e^+'^/') 

= [ (-d,{AUk){daUk + Ek dakm + AC/fc diaUk Vi " \[A\Jkf Va 

The symmetry of Uk implies: 

J^^ (^-d.iAUk) daUk + AUk diaUk vi - \{AVkf = 0. 



20 GILBERT WEINSTEIN AND LEI ZHANG 

In view of the equation ()3.ip , we also have 

/ d,{AUk)dS = 2hk{0) [ e^^K 

Substituting into (I4.13p . we obtain 

dahkiO) + 4hk{0)daM0) = 0{el^^) + 0(er'^')> a = 1,2, 3, 4. 

Now, if we choose r/2 + 1/2 < cj < 1, then 0(4^^) + 0{el'^^^) = 0(4^^), 
so that we can conclude: 

VhkiO) 



(4.14) 



hk{0) 



+ 4V(/.fc(0) 



4.2. The Vanishing Rate in g. In this subsection, we verify that (j4.14p 
leads to (jl.7p . For simplicity, we assume without loss of generality that 
the cut-off function x is supported in B{qik, 25) where 5 is small enough to 
guarantee that B{qijf_,26) are mutually disjoint. Indeed, this can be done 
since the left hand side of (jl.7p is invariant under any change of cut-off 
function x- Under this choice of cut-off, all the terms ViG{qik, qjk), j / h 
reduce to ViP{qik,qjk), so that it now suffices to show that: 

(4.15) 64^2 ^^p(^q^^^ ) _ 4VMQ^k) = -^M^ + 0(^/2). 

Indeed, 6"™''^''''=)/^ = 0(4^^) since \uk{qik) - Uk{qjk)\ < C for i / j, and 
furt hermore Vi(3{qik, qik) = '^2(3{qik,qik) since (3{x,y) = (3{y,x) by Lemma 
lA.ll in the Appendix. The remainder of this section is devoted to verify- 
ing (liT5]) . 

Taking the derivative with respect to ^ in (|2.5p and evaluating at qik, we 
have: 

(4.16) Vuk{qik) = 2[ ViH{qik,r])hk{7])e^'"''^^UVg{rj)+VMqik). 

J M 

where H{£^,'q) = — (l/Svr^) x(?') ^o?,dg^{£,,'r]). Similarly, for Uk we have: 

(4.17) Vuk{qik) = -^ [ V,{x{r)logdg{qik,v)phk{v)e^'"'^''^dVg{v) 



M 



+ / Vi/3(g,fc,r/)2/ifc(r?)e^"^('')dl4(r/)-V0fc((?ifc), 

J M 

where we used e^^* dVg^ = e'^'"= dVg. Let i?o(^,r7) = (l/87r2)x(r) \ogdg{i,ri), 
then we claim that: 

(4.18) \ViH{qik,r])-ViHQ{qik,ri)\ < Cdg{^,rj). 

Indeed, recall that w{qik) = and Vwk{qik) = 0. Thus, if fix ^ = q^k, and 
we let: 

fiv) = log dg(^,??) - log 4(^,77), 
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then Wkiri) = 0{\^ - r/p), and therefore |/(7?)| < C(|^ - It follows that 
V/(^) = 0. Now, by Appendix El Pgf{r]) is a bounded function, hence by 
elliptic theory, V'^f{r]) is bounded. We conclude that |V/(?7)| < C\S^ — rj\ 
from which (j4.18p follows. Next, since we have: 

it follows, by taking the difference of (|4.16p and (j4.17p . that: 
(4.19) 

= -VMQik)+ f Vip{qik,v)'2hk{v)e^''''^''^dVg{7]) -VMQik) + 0{ek) 
Jm 

Furthermore, we claim that: 
(4.20) 



N 

Vi[3{qik,v)2hk{ri)e*'"'^^UVg{r^) = IGtt^ V Vi/3(to, g,-^) + 0(el) 
Indeed, observe that: 



Vi/?((z,fc,77)2/ifc(r/)e4"'=(^)dF^ 
M 

N 



I V,P{q,k,v)2hk{v)e'''^'^''^dVg{rj) + 0{et) 



since, by (|1.6p . e *= = 0(e^,) on M \ l-)jLiB{qjk,6). In addition, for each 
j = l,..,N: 



Vi/3(to,r/)2/ifc(r/)e4"'=(^)dy3(7?) 



= Vi/?(g,fc,(?,fe) / 2hke^'"'dVg+ [ 0{\r, - qjk\)2hke^''HVg 

= 16^2^0(eD, 

where we used ()1.6p to estimate the first integral and a standard rescaling 
to estimate the second one. By combining ()4.14p with (I4.19P and ()4.20p . it 
follows that (I4.15|) holds. This completes the proof of the Theorem. 



Remark 4.1. Integrating (jl.2p . and using (|1.6p on the right hand side, one 
easily obtains: 

[ bkdVg = 87T^N + 0{el). 
Jm 

Appendix A. The Green's function for Pg 
Let G denote the Green's function of Pg as in (11. 3p : 

(A.l) f{i)-h=l G{^,v)Pgf{v)dVg{v), 
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where fg = Volg(M)~^ Jj^^ f dVg is the mean value of / with respect to g. 
Clearly, G is determined up to an arbitrary function of which we can fix 
by imposing the condition: 

(A.2) / G{^,v)dVg{7j) = 

J M 

The purpose of this appendix is to prove the following lemma, which is 
an improvement on a result of Chang- Yang [9j: 

Lemma A.l. Let {M,g) be a compact closed 4- dimensional manifold, and 
suppose Ker(Pg) = {constants}. Then the Green's function G{(,,r]) with 
respect to Pg can he written as: 

(A.3) = __i_^(^)iogr + /3(e,ry) 

where r = dg{^,r]) is the geodesic distance between ^ and rj, x(r) is a cut off 
function that is 1 on a neighborhood of ^ and vanishes outside B(^^, 6{^)/10j , 
and 6{S,) is the injectivity radius of {M,g) at ^. Furthermore, (3{(,,r]) G 
W'^''^{M X M), for any 1 < q < oo and satisfies: 

WPit -^w^-iiM) < C, uniformly in ^ £ M 
for some constant C = C{g,q). The principal part of G satisfies weakly: 

(A.4) P9,v[-^x{r)logdg{C,rj?j =5^ + E{tv) 

where E is a bounded function. Finally, we have G{S,,r]) = G{ri,^). 
Proof of Lemma \A . 1[ The weak form of (jA.ip is: 
(A.5) P,,,0(j,,,).5,__L_. 

Let g = e'^^g be a metric conformal to g such that in a neighborhood of ^ 
the normal coordinates with respect to g are conformal normal coordinates, 
i.e., det{g) = 1. Similar to (jA.Sp . we have for the Green's function Gi of Pg-. 

P-g,^Gl{i,ll) = 5^ ^ 



Volg(M)- 

By the conformal covariance of Pg this is equivalent to: 
Therefore since w{^) = 0, we have: 
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Note that the integral of the right hand side is zero, hence there exists 
F £ C^{M X M\D) such that: 

In particular, Gi = G + F, and the principal part of G and Gi differ by a 
bounded function. In the following, we focus on Gi. 

In conformal normal coordinates the point ^ corresponds to 0. We will 
identify y S T^M with its image under the exponential map. We first write 
Gi as in fO]) : 

(A.6) Gi(0,y) = -^x(^)logf + /5i(0,y), 

where r = \y\ = dg{0,y), and x is a cut-off function. We will show that the 
principal part H = — (l/87r^)x(f) log f of Gi satisfies weakly: 



(A.7) P-gH{0,y)=5o + Ei{0,y). 

where Ei is a bounded function. Since Pg/?i(0, y) = —Ei{0,y), this will 
imply by elliptic theory that /3i(0, •) G W^'i{M) C C3'"(M) for any 1 < 
q < oo. Here we assume that x = 1 in B{0, 6), x = in M \ B{S^, 26), and 
det(5) = 1 in B{0,25). 

Observe that since H{0, y) is radial and x is supported in a small neigh- 
borhood of where det(^) = 1, we have A?i/(0,y) = A'^H{0,y). Thus, it 
follows that for any smooth function (p: 

m = - f AlH{Q,y)ct>[y)dV-g(y)+ f HiO,y)Alcf>iy) dV-g{y). 

JM JM 

Clearly: 

/ A?if (0, y)^{y) dV-g (y) = [ A^i^ (0, y)0(y) dVg (y) 
JM JB2s\Bs 

and A?ff(0, y) is a bounded smooth function on B2S \ B^. Hence weakly, we 
have: 



(A.8) 



A?i7(0, y) = 6q + a, bounded function. 
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Also, {Pg — A'^)H{0, y) = Ax + a bounded function, where: 

= -^d^ [(^^m^m, - 2R^j{y)r'^ 9% 
= (^(^ldmR{y)Sm, - 29„(i?,,(y)rO) 9''^ 

lR{y)5^^-2Ri,{y)r'^d^-g'^^ 

= -^(^1 + ^2 + ^3). 

Here R, Rij, and Rijki are used to denote the scalar, Ricci, and Riemann 
curvatures of the metric g. The properties of conformal normal coordinates 
we used are listed in Appendix [Cj We estimate each of the three terms in 
the previous equation separately. 



Ai = (^^dmRiy)Sm, - 2i?i,-^(y)5™ - 2i?i,(y)9^r') 9 



ldmRiy)Smj - 2Rij4y) - 2Rij{y)dn,g"''^ ^ + 0(f) 

= -^R,,{0)yiyjr-^ + O{r). 

In the second equality above, we used g"^^ = 5mi+0(f^) and g''^ = 5ij+0{f'^). 
In the third equality, we used dmR{0) = and 2Rij^im{0) = Rjrn{0), as well 
as dmg"^^ = O(f^). It is easy to see that A2 = 0{f) since Rij{y) = 0(f), 
R{y) = O(r^), dmg"^^ = 0{f). Finally, we have 



^3 = {^-R{y)5^.j-2R,,{y)r'yg'^ 



-R{y)-g'^-2R,,{y)r''-g' 
^Rij{0)yiyjr^ - ^Riy)r'^ - 2Rij{y)6miSmjr'^ 



+ 4:Rij{y)6miSijyiymj r + 0{r) 
-^Rij{0)yiyjf^ + 2Rmi,ab{0)yaybyiymj + 0{f). 



ASYMPTOTIC EXPANSION FOR BUBBLING SOLUTIONS 25 

In the third equahty we used g'^^ = 5ij + O(f^), R{y) = 0{f'^), Rij{y) = 
0{f), and in the fourth equahty we used Rij^i{0)y^y^y^ = 0. Combining the 
estimates for Ai, A2 and A3, we obtain: 

A = ^ (^Rij{Q)yiyjf-'^ - 2Rij^im{0)yiyjyiy.mr~^^ + 0(f). 

Finally, this last estimate, together with (jA.SP yields (|A.7p as claimed earlier. 
Rewriting ()A.7p in integral form, we get: 

(A.9) m= [ H{0,y)P-g<P{y)dVg{y)- [ Ei{0,y)<P{y) dVg{y). 
Jm J m 

Note that i?i(0,-) is supported in B{S,,25) \ In this small neigh- 

borhood, let r] = exp^(y), where exp^ : T^{M) — > M is the exponential map 
with respect to the metric g. Thus, recalling that the principal parts of G 
and Gi differ only by a bounded function, we may rewrite (|A.9h as: 

m = - [ AxW logdgiC,v)P9mdVg{rj) 
(A.IO) ■^''^^'''^ 

JB{^,25)\Bi(,S) 

Written weakly, this reads: 

Pg ^~2''^(^) ^ogdg{S^,r])^ = + a bounded function, 

as claimed in Lemma lA.ll 

Finally we show r]) = G{i], ^). This part of the proof is similar to the 
proof in [3l page 108]. We include it here for completeness. We first claim: 

(A. 11) / G{ri,(,)dVg{i]) = constant. 

Indeed, let /(O = Jm G{ri, dVgir]), then for any (j) G C7^(M), 

/ f{0P9,imdVg{.0= I ([ Girj,C)Pg,^mdVgiC)) dVgirj) 

JM JM \JM J 

= [ {m-^g)dVg{r,)=Q. 
JM 

Hence Pgf = in weakly, and since Ker(Pp) = {constants}, we have proved 
([AlT]) . Next for any (() e C^(M), 

m -4>g= f G{rj, OPa,^m dVgiO. 
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hence for any ip E C^(M): 

/ Pg,^mmdVg{0= [ mP9,v^{ri)dVg{r^) 
JM JM 

JM JM 

JM JM 

Since Ker(Pg) = {constants}, we obtain 

H0= / G{'n,^)Pg^ni'{v)dVg{r]) + constant. 

JM 

Thus, using also the definition of G(^,r/), we obtain 

I (G(e,r/) - Gi7],0)Pg,^Hv) dVg{7]) = constant. 

JM 

Integrating with respect to and using ()A.2p and (|A.lip . we see that the 
constant on the right hand side above is 0. Thus, again by (jPj) we have 

G{tr])-Gi7j,0 = G. 

where G is independent of rj. Since the left-hand side is clearly antisymmetric 
with respect to ^ and ry, we obtain G{S,,rj) = G{rj,S). This completes the 
proof of Lemma lA.li □ 

Appendix B. Comparison between d^{y,z) and \y — z\ 

In this subsection we establish the following estimate: 
(B.l) 

Di{^og\y - z\ - log dg^^{y,z)) = 0(e||yp"^), for \z\ < -y and \y\ < Sis^'^, 

for j = 1,2,3. We recall that gk = £^'^4>*9k is the blow-up metric, and 
we identify y,z €z TqM with expy and expz respectively, where exp is 
the exponential map at the origin with respect to the metric ^fc. Thus, 
dgk{y,z) = d§^.(expy,expz). 
Fix y and consider 

2 

f{z) =log|y-2:| - \ogdg^{y,z), on \z\ < -\y\, 

as a function of z. We shall obtain the following estimate: 

(B.2) \Dif{z)\<Cel\y\^-^, for |z| < |y|/2, j = 1, 2, 3. 

Then since L>^/(z) = (-1)J\d|/(z), (IBH) follows. 

Let R, Rij, and Rijki respectively denote the scalar, Ricci, and Riemann 
curvature of g, where we suppress the subscript k. By the definitions of g 
and FCji^iz) one obtains easily that 

v^i?;.^,(z) = o(4+^), j = o,i,2. 
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As a consequence, using R{0) = |V-R(0)| = Rij{0) = 0, we have further: 

(B.3) R{z) = 0{st\z\^), R,,{z) = Oiel\z\). 

We also note the following simple estimate on df. 

(B.4) \Di{logdg{y,z))\<C\y-z\-^, j = l,2,3,4. 

We shall now derive an estimate on A?/(z). By ()A.4p and the definition 
of ^, we have: 

(B.5) Pg,Aogd^g{y,z) = 0{et). \z\<'^\y\, \y\ < Sie^\ 

Next, using (|B.3P and ()B.4p . we can estimate the term: 
(Pg-A|)logdg(y,z) 

= dm (5™ (^lR{z)g^J - 2R,,{z)^ g^^dj{logd^g{y,z))^ = 0{el\y\~'). 

Combining this with (|B.5p . we get: 

(B.6) A|,(logdg(y,z)) = 0{el\y\-'), \z\ < ^\y\, \y\ < 5,e^\ 

Finally, we consider the term A? ^(log \y — z\). Since A^(log \y — z\) = 0, 
it suffices to estimate A? ^ — A^. For any function u, we have, by direct 
computation: 

(B.7) Al,u = r'pd.jabu + 2dijau{d^r''9'' + r'dbP) 

where we used det(^) = 1. Using the expansion of gij{z): 

= Sij + yiRpijg{o)z'^ + 0(4|z|3). 
where z*-' = z^z^ , and replacing u by log \y — z\ in (|B.7p . we obtain: 
A|,(log \y - z\) = 0{el\y\-'), \z\ <^\y\, \y\ < 6^8^^ 
and consequently 

(B.8) ^J{z) = 0{el\y\-^), \z\ <\\y\, |y| < ^ 

An estimate on the L°° norm of f{z) is easily obtained: 

(B.9) d5(y, z) = j' ^g,,{t)x[{t)x'^{t)dt = \y-z\{l + 0{el{\y\^ + \z\^))). 

Note that here we didn't need the assumption: \z\ < ||?/|. From ()B.9p . we 
clearly have 

(B.IO) /(z)=0(4|y|2) \z\<^\y\, \y\<6^e^\ 
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With (jB.Sp and (1B.10|1 . we can apply the standard rescahng argument and 
elhptic theory to derive ()B.2|) . 

Appendix C. Conformal normal coordinates 

In this subsection we hst some well known facts for convenience. Let g be 
a metric on M and let p S M be a point. We emphasize that g and u are 
unrelated to the corresponding quantities in other sections. Suppose that 
in normal coordinates around p the metric g satisfies det{g) = 1, i.e., those 
normal coordinates are conformal normal coordinates at p. We will denote 
p as and consider the properties of g( in a neighborhood of 0. 

First, since det{g) = 1, the Laplacian is given by: 



(C.l) AgU = djig'^dju) = djg'^djU + g'^diju. 

The second term PgU — AgU of the Paneitz operator is given by: 

(C.2) divg - 2Ric^ du^=dm - '^Rij^ diug^^ g 

Further properties of conformal normal coordinates include: 

(C.3) i?y(0)=0 

(C.4) i?,,- fc(0) + R,kM + i?fe,,,(0) = 

(C.5) Vi?(0) = 

(C.6) AR{Q) = - V(0)|' 

D 

2 

(C.7) R(ij^ki){{)) + -R(pijmi^)Rpklm){0) = 

(C.8) Rpijg^p{0) = Riqj{0) — Rij^q{0). 



Here, we use the round brackets as the customary notation for the symmetric 
part with respect to all indices not contracted within a certain scope. 
In normal (not necessarily conformal normal) coordinates, there holds: 

(C.9) g^biO = Sab + iRaijbm'' + lRa^Jb,me'^ + 0{r^), T = \C\ 

and consequently, the same expansion holds for the inverse: 

(C.IO) 5"'(e) = Sab - lRa^Jbme' - lRaijb,km''' + Oir*). 

Here ^^-^ denotes C C'^ i etc. 
Taking a derivative: 

(C.ll) d,g''\0 = -lRa(c^)biO)e-l{2Raic^)bA^) + Ra^,b,M)e' + Oir^). 

Contracting over to a and c and using (jC.Sh : 
(C.12) dag'^'iO = - J(2i?.6,,(0) - R^,,b{0))e' + 0{r^). 
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Taking another derivative in (IC.llll : 

(C.13) d^dQ^Hi) = -\Ra{cd)m 

Contracting over a and c and using (1C.8|1 and (|C.3|) : 
(C.14) dadg^^i) = \R^d,me + 0{t^). 

Appendix D. A Pohozaev Identity 
In this appendix, we derive a Pohozaev identity for the equation 
(D.l) Pgu + 2b = 

Throughout, we assume that det{g) = 1 over 0, which we take to be a bah 
centered at 0. 

First, multiplying the right hand side of ()D.ip by S,^diU and integrating 
by parts, we have: 



/ 



^ Jan Jn ^ 



where Vi is the unit normal to the boundary dVl. Also note that we omitted 
to write here, as we will throughout this appendix, the volume element dVg 
in the integral over O, and the area element dAg in the integral over dO.. 

Next, we consider the term A^it on the right hand side of (jP.ip . Again, 
multiplying by the same factor ^ • Vu and integrating by parts, we get: 

/ A2u(^Vn)= / di{g'^d,{^gu))edku 
Jn Jn 

= [ g'^djiAgu)C''dkUUi- [ g'^dj{Agu)d,{edku) 
Jan Jn 

= h- I g''d,{Agu)d^u- [ g'^dj{Agu)tdikU 
Jn Jn 

= h- g'^AgudiUv,j+ {Aguf- g'^AgU^^dikUVj 
Jan Jn Jan 

+ [ Agudjig'^^diku) 
Jn 

= h-h + I {^guf -h+ [ g''Agud,,u + [ Agutdjig'^d^ku). 
Jn Jn Jn 

Here /i, I2 and Is are boundary integrals. In order to compute the last two 
terms we now note: 

g'^diju = AgU - dig'^dju 
djig'^diku) = dk{Agu) - dikg'^dju - dkg'^diju. 
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These imply: 

/ g'^AgudijU= [ (Aguf- [ Agudig'^dju= [ {Aguf - B3, 
Jn Jci Jci Jci 

and 



AguCdjig'^diku) 

JQ Jo, Jn 

= \j (^ u){Aguf - 2 / {Aguf -B,- B2. 
^ Jan Jn 

Thus, we conclude: 

f Alu{^ ■ Vn) =Ii-l2-h + hi -B1-B2- B3. 

Finally, we consider the second term in the Paneitz operator. As before, 
we multiply by ^^d^u and integrate by parts: 

= Ci - ^ (^^^^^ ~ 2i?i, (o) g'^diu {dmu + 

Since Rij{0) = 0, i?(0) = 0, and diR{0) = 0, we see that Ci above can be 
written as 

Ci= f i-2Rij,i{0) djudku^ e + 0{t'')\DuW T=\i\. 

JdQ 

We also estimate the second term: 

J^g""' (^lRi09ij-2Rij{0^ diug'^{dmu + edmku) 

= I {{0{r^)-2Ri,^s{^)i')^lug^\^iuH''^^ku)+0{r'^) \Du\ {\Du\+r\D\\)) 
Jn 

= I ((0(r2) - 2i?i,- ,(0) r ) dju {d,u+e diku) + 0{r^) \Du\ {\Du\+r\D''u\)) 
Jn 

' {-2Rij,m ^^dju{diu + diku) + 0{r^)\Duf + 0{r^)\D^u\) . 
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Finally, putting all the estimates together, we obtain the final form of our 
Pohozaev Identity: 

(D.2) J + a^/ie^"^ 

Jan 

+ g'^ Agu e diku \ (Aguf C ^^i) 

^l[Agu0.g^^0,u^eAgua..g^^a,uU^Agu8.g^^8,u-2K0.u) 

+ 2 / (Rij,i{0) djudkuC' i^i + 0{r^)\Du\A 
Jan ^ ' 

- f {2Rij,ii0){djudiue + djudikue^') + O{r^)\Du\^ + O{r^)\D^u\). 
Jn 
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